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Abstract. - We study a system of two symmetrical capacitively coupled quantum dots, 
each coupled to its own metallic lead, focusing on its evolution as a function of the gate 
voltage applied to each dot. Using the numerical renormalization group and poor man's scaling 
techniques, the low-energy Kondo scale of the model is shown to vary significantly with the 
gate voltage, being exponentially small when spin and pseudospin degrees of freedom dominate; 
but increasing to much larger values when the gate voltage is tuned close to the edges of the 
Coulomb blockade staircase where low-energy charge-fluctuations also enter, leading thereby 
to correlated electron physics on energy /temperature scales more accessible to experiment. 
This range of behaviour is also shown to be manifest strongly in single-particle dynamics and 
electron transport through each dot. 



Introduction. - Recent years have seen intense investigation of electron transport through 
semiconducting quantum dot devices [1] , leading in turn to a strong resurgence of interest in 
basic Kondo physics [2]. The classic example — namely the spin- 1/2 Kondo effect [3] in which 
the magnetic moment of an odd-electron dot is quenched by coupling to metallic leads — was 
observed in quantum dots almost a decade ago. Since then, a continuing goal has been to 
understand how coupled, multiple quantum dot systems may lead to variants of the Kondo 
effect involving both spin and orbital degrees of freedom. 

In this Letter we consider probably the simplest example of such — a symmetrical double 
quantum dot system in which the interdot coupling is capacitive. Experimental realisations 
of such devices have appeared in the literature [4-6] , and various aspects of the rich inherent 
physics have been uncovered in a number of theoretical papers, see e.g. [7-13]. Here we 
elucidate theoretically a key underlying issue: the effect on transport of sweeping the dot 
energy levels through a wide range of values by means of a suitably applied gate voltage. 

The low-temperature transport through each dot is heavily dependent on the magnitude of 
the gate voltage: both the characteristic low-energy Kondo scale of the system, and the zero- 
bias differential conductance, vary significantly as the energy levels of the dots are lowered 
from the Fermi level down. If the gate voltage is tuned in such a way that only low-energy 
spin, or spin/orbital-pseudospin, degrees of freedom are relevant, then the Kondo scale is 
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found (as one usually expects) to be exponentially-small. By contrast, at points where low- 
energy charge fluctuations can arise, we show that the Kondo scale is much larger — even when 
the system is intrinsically strongly correlated; leading thereby to correlated-electron physics 
on energy scales more amenable to experimental interrogation. We consider, and analyse, this 
behaviour within a renormalization group (RG) framework. 

Model and atomic limit. - The capacitively-coupled double dot model we consider con- 
sists of two correlated Anderson impurities coupled by an interdot Coulomb repulsion, with 
each dot i — L, R coupled to its own separate lead; viz 

H^Y. '^koL^ak,^ + ^("Lo-c^^ + li-c-) + Yl^^^' + Uh,^n,i) + U'tilur (1) 

with flier = cl^Cia- and fii — flier- The first term here describes the two metallic leads 
(with i = L,R again). The second term represents the tunnel coupling of each lead to its 
corresponding dot, and the final terms define the energies of the double dot itself (including 
both intra- and interdot Coulomb repulsions U and U'). The beauty of a typical experimental 
quantum dot device is that the parameters entering its effective Hamiltonian are readily 
adjusted. In particular the dot levels e can be tuned over a wide energy range by varying the 
gate voltage. It is this on which we focus here, considering eq. ([J) as a function of e when 
U' < U. We also consider e < only, i.e. e — — |e|, since the alternative corresponds to a 
'frozen impurity' regime with the dots in essence unoccupied and the basic physics relatively 
straightforward. Moreover, the behaviour of eq. ([l]) for |e| above the particle- hole symmetric 
point |e| = U' + U/2 is obtainable from that below it by a simple particle-hole transformation; 
we can thus restrict our analysis to the range < \e\ < U' + U/2. 

In the V = ('atomic') limit the dots are decoupled from their leads. The total dot 
occupation number n versus |ej then follows the familiar Coulomb blockade (CB) 'staircase' 
pattern, increasing stepwise from n = to n = 4 as |e| is increased. The parameter range 
of interest involves only the n = 0, 1 and 2 steps, the edges between them lying at |e| = 
(n = <-> 1) and |e| — U' {n — 1 ^ 2). For all points along the n — 1 step, the degenerate 
(rii, nji) = (0, 1) and (1, 0) configurations lie lowest in energy, whereas on the n = 2 step the 
ground state is (1, 1); and at the n = 1 ^ 2 CB edge these states are all degenerate. 

In what follows we elucidate the behaviour of eq. ([T]) with increasing |e| when the interacting 
dots are connected to their leads. We show that for all U' < U, the 'atomic limit' degrees 
of freedom described above are always quenched as T ^ 0, such that the ground state is a 



non-degenerate Fermi liquid phase throughoui ll'^')! More central, however, is how the nature 



of this quenching depends on the gate voltage |e|. Near the midpoints of the CB steps, 
the dot degrees of freedom (be they spin or a combination of spin and orbital-pseudospin) 
are Kondo-quenched on exponentially small temperature scales. Near the edges of the steps 
however, where charge fluctuations are important, the behaviour is markedly different. The 
Kondo temperature Tk is found to be much larger, on the order of the dot-lead hybridisation 
parameter F = irpV^ (with p the host /lead density of states) . Since Tk sets the natural scale 
below which the many-body physics arises, its enhancement means that non-trivial coherent 
electron transport occurs over a much larger temperature window. 

NRG and fixed points. - Results discussed in the next section have been obtained using 
Wilson's powerful numerical renormalization group (NRG) approach [14, 15], technical discus- 
sion of how to apply which to the Hamiltonian eq. ([T]) can be found in [12]. NRG is of course 

(^)It is known [11, 12] that in the n = 2 regime for a critical U[. > U, the model undergoes a quantum phase 
transition to a non-Fermi-liquid, doubly-degenerate 'charge ordered' phase. 
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Table I - NRG fixed points discussed in this work. For detailed discussion, see text. 
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00 











N/A (since F = 00) 
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an iterative procedure that converges ultimately to a stable fixed point [14, 15], the nature of 
which indicates the structure of the ground state; and the NRG flows that lead to it (via the 
unstable fixed points) can be used to obtain static and dynamic properties on all important 
temperature and energy scales respectively. 

Our results will be interpreted in terms of six distinct fixed points, the Hamiltonians for 
which may be obtained simply by setting each of the bare parameters of eq. ([T]) to either or 
00. The simplest is the unstable free-orbital (FO) fixed point [14], where e = U — U' ^ T — 0. 
Four more unstable fixed points arise from F = and U ^ 00 (with fixed ratios oi e/U and 
U' /U), each corresponding to one of the atomic limit states discussed above. Finally there is 
the stable strong-coupling (SC) fixed point in which F = 00, to which all NRG flows ultimately 
tend. Table [J lists all these fixed points, the conditions under which they obtain, the allowed 
dot configurations and their degeneracies. 

As seen in table U the unstable fixed points divide into two sets. The first comprises 
the 'local moment' (LM) fixed points, in which the total dot charge is fixed and hence the 
only possible dot degrees-of-freedom are those of spin and orbital pseudospin. The second 
set contains the 'valence fluctuation' (VF) fixed points, in which charge fluctuations of the 
double dot also enter, ft is the existence of these two distinct types of fixed point that leads 
to interesting physics, as now explained. With decreasing temperature/energy scale, all RG 
flows begin close to the FO fixed point (we take the host bandwidth D to be the largest energy 
scale of the model). But for fixed U and U' , the subsequent flows then depend solely on the 
magnitude of [e[ (as controlled by a gate voltage); such that the flow from the final unstable 
fixed point to the stable SC fixed point determines the characteristic low-energy scale. 

For most values of [e[ — sufficiently far from the CB edges — the NRG fiows are found to 
pass essentially from FO to one of the LM fixed points, before crossing over to the stable 
SC fixed point where the dot local moments are quenched. The effective low-energy models 
here are thus Kondoesque and the low-energy Kondo scale Tr- on which the SC fixed point is 
reached is exponentially-small in U/T and U' /T (as considered further below). 

If however [e[ is instead tuned to within ©(F) of a CB edge, the VF fixed points also come 
into play and the resulting physics is quite different. This is most easily seen at one of the CB 
edges itself, where there is no sign of any LM fixed points at all. The RG flow instead goes 
from FO to SC via VF, and the Kondo scale is found to be a simple multiple of F, virtually 
independent of the bare U and U' . So by tuning the system from one CB step to another, one 
can produce dramatic changes in its Kondo scale and low-energy degrees-of-freedom, which 
are of course evident in both statics and dynamics. 

Before discussing numerical results, we make one further remark. The Kondo physics — 
and in particular the exponentially small Tk — associated with the LM — > SC crossover is 
a well-known manifestation of strong electron correlation. But while the low-energy scales 
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associated with the CB edges depend only on the bare energy scale F, the corresponding 
electron dynamics are nonetheless still correlated. Indeed for U /T sufficiently large, the only 
region in which electron correlation is not important is e > F. As we shall see below, this 
means that the physics associated with crossing the n — \ ^ 2 CB edge is distinct from that 
arising near 0^1 [7]: only the former involves correlated physics throughout, and it is this 
to which we direct the majority of our attention. 

Results. - We first present some static properties, beginning with the entropy. Since 
each of the fixed points in table [J has a well defined dot degeneracy, the contribution to the 
entropy of the system from the dots alone, denoted by 5'imp, can be used to determine how 
the NRG flows depend on the bare parameters. In what follows, it is convenient to work with 
the dimensionless parameters |e| = |e|/7rF, U = U/ttT and U' — C/'/ttF. 

Figure [TJa) shows 5'imp as a function of temperature for fixed [U ,U') — (14,12) at four 
representative points along the |e| line. These are the three points of 'symmetry' at |e| — U' /2 
(centre of n = 1 CB step, solid line), \e\ = U' {n — 1 ^ 2 CB edge, long dashes) and 
|e| = [/' + U /2 (centre of n = 2 CB step, dot-dashed); plus a fourth point at |e| = 11 (short 
dashes) which illustrates the behaviour slightly away from the n = 1 ^ 2 CB edge. In all 
cases, the high-temperature In 16 free-orbital behaviour is seen to cross over ultimately to 
the singlet strong-coupling regime below some characteristic low temperature scale, but at 
intermediate temperatures the physics is clearly different in each case. 

The |e| = U' /2 and U' + U/2 curves, corresponding to the centres of the n = 1 and n = 2 
CB steps respectively, show the generic Kondoesque behaviour described above. The long 
intermediate plateaus of In 4 each indicate a flow that passes very close to one of the two 
LM fixed points in table [J Specifically, for |e| = t7'/2 the fixed point is LMf^^^-* and the 
low-energy effective model is the one-electron SU{A) spin/orbital Kondo model [8,9]. For 
|e| — [7' + {//2 by contrast, it is the LM^'^^^'^ fixed point that is approached, and at low 
energies one observes an effective pair of uncoupled SU{2) spin-Kondo models [11]. 

The |e] = U' entropy curve in fig. [Ija) illustrates the distinct behaviour that arises when 
the RG flow passes close to a VF fixed point — in this case the VF12 fixed point of table [B 
No local moment plateaus are observed; low-energy charge fluctuations are now possible and 
the entropy instead drops to zero on a much higher temperature scale characterised by the 
hybridisation parameter F. When one is sufficiently close to, but not at, the CB edge at 
|e| = [/', both VF and LM fixed points are observed, as in the ]e| = 11 (= U' — 1) curve in 
fig- Ela); here the entropy follows the |e] = tj' form at relatively high temperatures before 
crossing over to a characteristic In 4 LM shoulder and then finally zero as T — ^ 0. 

To illustrate charge fluctuations near the 1^2 CB edge, fig. [TJb) shows the T = charge 
susceptibility Dx^imp(O) 1^1 (with 4Tx+(T) ^ {[N - {N)f) and N the total charge operator 
[12]). A significant peak occurs at the CB edge at \i\ = U' = 12, refiecting the availability here 
of the low-energy charge degrees-of- freedom. Away from this point however, the low-energy 
effective model involves only spin or spin/orbital degrees of freedom and xi^imp(O) very 
small (vanishing as t7 ^ 00). The HWHM of the peak gives an estimate of the range over 
which low-energy charge fluctuations are significant: here we see it is 0(1), and hence the 
charge degrees-of-freedom are important only when |e] is within C(F) of the CB edge. 

We now analyse the Kondo temperature Tk as a function of ]e] (with Tr- taken to be the 
temperature at which the impurity spin susceptibility satisfies kBTxs,\m-p{T) j {gfiBY — 0.1). 
As seen already in fig. [DJa), Tk varies by many orders of magnitude as e is lowered: this 
is shown more clearly in fig. [2ja) where we plot IuTk/D versus —c/ttT for the same fixed 
(f/, f/') = (14, 12) (solid line). Also shown is the total dot charge n (dotted line), showing the 
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characteristic 'rounding' of the CB staircase on coupUng to the leads. It is quite clear that 
Tk is a minimum when n ~ 1 or 2, and a maximum of 0(T) at the 1 <-> 2 CB edge when 
n can fluctuate easily between the two. Note that Tk also increases as one moves from the 
n ~ 1 SU(4:) Kondo regime towards the -s-^ 1 CB edge; but in contrast to the behaviour 
just described continues to increase as e becomes positive, reflecting the essential absence of 
electron correlation in the n ~ regime and a Kondo scale naturally on the order of F. 

The dependence of Tk on |e| can be understood more precisely from the effective Kondo 
models in the n = 1 and n = 2 regimes, obtained by straightforward Schrieffer- Wolff transfor- 
mations and then analysed using Anderson's "poor man's" scaling [3]. We find in both cases 
that the leading exponential form of Tk/ D is exp[— (2/9j)^^], with pj given by 

, {2/\~e\-{\~e\-Ur'~m~U)-' ior \~e\ U' /2 

\{\A-U'y^ -{\i\~U' -il)-^ for |e| - [/' + C//2 ^' 

in the strongly correlated {U 3> F) limit. These results for Tk are plotted as dashed lines in 
fig- m^), for the same fixed U and U' , each scaled by a single multiplicative factor onto the 
NRG result at the appropriate minimum. The strong agreement with numerical and analytical 
results over a wide range of |e|, provides further confirmation that the low-energy physics of 
the dots is essentially captured by Kondo models of fixed n, except when |e| is within 0{r) 
of a CB edges and charge fluctuations become important. 

In Figure [2jb) we also plot IuTk/D versus |e| for a range of different coupling strengths. 
The key point here is that while the scales in the Kondo regimes clearly vary significantly 
with all the bare parameters of the model (as seen from eq. ([2])), those associated with the VF 
points are largely dependent on the hybridisation strength F only. For this reason, one would 
expect all the low-temperature thermodynamics of the model to be essentially independent of 
U and U' near the CB edges, and we have indeed found this to be true. 

We now turn to dynamics, focusing on the local single-particle spectrum Di{uj) defined in 
terms of the retarded Green function Gi{t) = —i6{t){{c^^{t), cl^}) via Di{uj) = — ImG'i(w), 
since this quantity probes the charge fluctuations that control transport across a given dot i 
(the model considered, where each dot hybridises to its respective lead, is formally equivalent 
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Fig. 2 - Variation of Tk with [^I. (a) Tk/T vs \i\ for {U, U') = (14, 12): solid line from NRG, dashed 
lines from poor man's scaling. The dotted line shows the average total dot occupation number 
n (right-hand scale), (b) NRG Tk versus |e| for different inter- and intradot coupling strengths: 
{ty,U') = (10,8) (solid line), (12, 10) (long dashed), (14, 12) (short dashed) and (18, 12) (dot dashed). 
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to a 4-lead device where each dot hybridises to two leads with strength r/2). As shown 
in [12, 16], the T = zero-bias differential conductance (G) of each dot is related exactly 
to the zero-frequency spectrum by G/(2e^//i) = irVDiifi) (and the spectrum at non-zero lo 
provides an approximation to the conductance at source-drain bias V = uj/e). 

The NRG technique can be used in practice to calculate Di{uj) on all frequency scales of 
interest [17]. Near the centre of each CB step, we find that the high-frequency spectrum is 
dominated by its Hubbard satellites, corresponding physically to sequential tunneUing through 
the dot. Since these features are readily explained simply as broadened versions of the poles 
arising in the atomic hmit = 0, we do not discuss them further here. More interesting is 
the low-frequency form of the spectrum, since this corresponds to coherent tunneUing through 
the dot mediated by the Kondo effect. In figures [2la),(b), we illustrate the evolution of the 
low-energy Kondo resonance with increasing |e| for fixed {U,U') = (10,8). Figure[3]Ja) shows 
the spectrum for |e| = 4, 6, 7 and 8 (dot-dashed, dotted, dashed and solid respectively), 
and fig. mb) shows |e| = 10, 11 and 13 (dashed, dotted and solid). Note the very different 
frequency scales in each case. The width of the resonance, being proportional to Tr-, is seen 
to first increase with |e| as one moves from the centre of the rt = 1 CB step (|e| — U' /2 = 4) 
to the maximum at the 1 <-> 2 CB edge |e| = U' = 8, where Tk ^ 0{T). Thereafter, the 
resonance narrows until the particle-hole symmetric point at |e| = U' + U /2 = 13 where Tk 
is a minimum. In addition to the changing width of the resonance with |e|, there is also a 
definite increase in asymmetry near the 1^2 CB edge, since here the addition and removal 
of low-energy electrons to and from the ground state clearly occur with different weights. 

Finally, the zero- frequency spectral density T:TDi{0), i.e. the zero-bias differential con- 
ductance, is shown as a function of |e| in fig. [3jc). Its stepwise increase is readily explained 
by Fermi-liquid theory, as one can show that 7rTDi{0) = sin^(7rni/2) with rii (= n/2) the 
occupation number of dot i, whence T:TDi{0) itself mimics the characteristic CB staircase 
form. For T <IC Tk, the zero-bias differential conductance provides as such a direct means of 
measuring the average occupation number of each dot. 

Conclusion. - The capacitively coupled double dot system described by eq. ([T]) shows a 
rich range of behaviour as the energy levels of the dots are varied, which we have elucidated by 
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Fig. 3 - Spectral density of dot i for fixed {U, U') = (10, 8). (a) and (b) show the evolution of TTrDi{u}) 
with increasing |e|: specifically, |e| = 4, 6 and 7 and 8 (= U') [(a): dot-dashed, dotted, dashed and 
solid respectively]; and 10, 11 and 13 [(b): dashed, dotted and solid], (c) T = zero-bias conductance 
G/{2e'^/h) = 7rrA(0), as a function of 




a combination of NRG and poor man's scaling techniques. By applying a suitable gate voltage, 
the system can be tuned to display either an SU{2) x SU{2) spin-Kondo effect {n = 2), an 
SU (4) spin/orbital Kendo effect {n = 1), or correlated mixed- valence physics where the Kondo 
scale is substantially enhanced and non-trivial coherent electron transport thus arises over a 
much wider range of temperatures. It can moreover be shown that this essential underlying 
physics is robust to inclusion of a direct interdot hopping ('t'), provided \t\ < Tk- 

* * * 
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